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II.  AlUTf  #  C  1 

Equations  of  stata,  or  constitutive  aquations,  provide  a  complete  mechanical  as 
wall  as  thermodynamic  description  of  material.  An  aquation  of  state  relates  stress  or 
pressure,  strain  or  density,  and  one  or  more  thermodynamic  variables,  such  as  Internal 
energy,  or  temperature,  for  s  specific  materiel.  ^Equations  of  state  for  earth  materials 
ore  necessary  for  the  calculation  of  Intense  stress*weve  propagation  and  seismic  coupling 
|  resulting  from  explosions  or  Impacts  on  or  within  the  earth's  crust,  Beacuse  of  the  wide 
range  of  dynamic  pressures  which  ere  Induced  within  the  vicinity  of  a  nuclear  explosion, 
the  equations  of  state  which  ere  of  Interest  will  describe  the  response  of  a  medium  over 
a  range  of  conditions  from  thousands  of  kllobars  and  tons  of  thousands  of  dsgrees  down 
to  s  fraction  of  e  bar  end  ambient  temperature;  the  letter  conditions  correspond  to 
large  distances  from  the  disturbance. 


The  present  program  has  been  cone a reed  with  the  ansi ye Is  of  e  large  body  of 
published  data  obtained  hy  shock-wave  techniques  which  describe  the  dynamic  behavior  of 
e  npmber  of  earth  materials  below  2  megabers.  Most  of  those  data  are  expressed  In  tame 
of  Hugonlot  curves,  hcvever  the  e»*hors  have  used  additional  data  such  as  composition 
of  the  rocks  or  minerals  lnuol  i  he  nondynamic  properties*  and  attempted  to  examine 
the  systematise  of  the  equation*  *  *  ste  of  these  materials.  The  authors  have  also 
reduced  the  shock-wave  data  which,  moral,  are  In  terms  of  a  pressure- vo  1mm- 
energy  relation  to  the  mere  convent**  pressure- vo  1«m- temperature  variables. 
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I.  Introduction 

Equations  of  state,  or  constitutive  equations,  provide  a  complete 
mechanical  as  well  as  thermodynamic  description  of  material.  An  equation 
of  state  relates  stress  or  pressure,  strain  or  density,  and  one  or  more 
thermodynamic  variables,  such  as  internal  energy,  or  temperature,  for 
a  specific  material.  Equations  of  state  for  earth  materials  are 
necessary  for  the  calculation  of  intense  stress-wave  propagation  and 
seismic  coupling  resulting  from  explosions  or  impacts  on  or  within  the 
earth's  crust.  Because  of  the  wide  range  of  dynamic  pressures  which 
are  induced  within  the  vicinity  of  a  nuclear  explosion,  the  equations 
of  state  which  are  of  interest  will  describe  the  response  of  a 
medium  over  a  range  of  conditions  from  thousands  of  kilobars  and  tens 
of  thousands  of  degrees  down  to  a  fraction  of  a  bar  and  ambient 
temperature;  the  latter  conditions  correspond  to  large  distances  from 
the  disturbance. 

The  present  program  has  been  concerned  with  the  analysis  of  a 
large  body  of  published  data  obtained  by  shock-wave  techniques  which 
describe  the  dynamic  behavior  of  a  number  of  earth  materials  below 
2  megabars.  Since  most  of  these  data  are  expressed  in  terms  of  Hugoniot 
curves^  we  have  used  additional  data  such  as  composition  of  the  rocks 
or  minerals  involved,  thermodynamic  properties,  elastic  properties, 
crystal  structure  and  strength  measurements,  and  attempted  to  examine 
•  the  systematics  of  the  equations  of  state  of  these  materials.  We 

have  also  reduced  the  shock-wave  data  which,  in  general,  are  in  terms 


of  a  pressure-volume-energy  relation  to  the  more  conventional  pressure- 
volume-temperature  variables. 

II.  A  Review  of  Equation  of  State  Data  for  Earth  Materials  Pertinent 
to  Seismic  oupling 

Equation- of- state  data  are  obtained  in  the  very  high  pressure 
range  —  the  hydrodynamic  regime  —  with  both  explosive  and  hypervelocity 
Impact  shock  wave  techniques.  Because  the  magnitude  of  the  principal 
stresses  (for  example,  the  radial  stress  In  spherical  shock  propagation) 
far  exceeds  any  possible  stress  differences  supportable  by  the  medium, 
a  fluid- like  behavior  is  usually  assumed  in  this  regime.  Although  the 
formulation  of  the  equation  of  state  in  this  regime  lacks  the 
complications  arising  from  material  strength,  the  fact  that  many 
minerals,  In  particular  the  silicates  (and  rocks  containing  silicate 
minerals)  undergo  solid-solid  (polymorphic)  transitions  (Boyd,  1964 
and  Rlngwood,  1969)  turns  the  formulation  of  a  complete  equation  of 
state  Into  a  challenging  problem.  In  addition,  as  a  result  of  the 
high  temperatures  which  accompany  the  high  pressures  In  the  hydrodynamic 
experiments,  descriptions  of  partial  or  complete  transformation  of 
rocks  or  minerals  into  liquids  and/or  gases  are  required. 

At  distances  far  enough  from  the  source  that  the  dynamic  stress 
levels  are  below  the  dynamic  yield  strength  of  the  medium,  an  elastic 
type  of  behavior  Is  usually  assumed,  as  a  first  approximation,  for 
the  material.  In  the  elastic  regime  the  equations  of  state  are  usually 
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derived  from  laboratory  low  amplitude  shock,  acoustic  (ultrasonic), 
and  quasi-static  pressure-volume  experimental  data. 

Over  a  range  of  intermediate  stress  levels  (e.g.  tens  to  hundreds 
of  kllobars) ,  dynamic  yielding  of  rocks  and  stlnerals  la  uaually 
encountered.  Experimentally  it  is  observed  that  for  a  range  of  ahock 
stresses,  a  single  shock  wave  transmitted  into  the  material  will  break 
into  an  elastic  and  a  deformatlonal  shock  front.  These  shock  fronts 
are  thought  to  represent  the  onset  of  dynamic  uniaxial  compresalons ; 
these  are  achieved  by  essentially  elastic  and  deformatlonal  strain, 
respectively.  Of  special  Importance  In  this  regime  is  the  relation 
of  stresses  parallel  and  perpendicular  to  the  stress  wave  propagation 
direction  during  and  after  the  material  has  been  encompassed  by  the 
deformatlonal  shock  front. 

Relevant  to  the  development  of  complete  equations  of  state  for 
rocks  and  minerals  are  contributions  from  the  following  sources: 

1.  Shock  wave  experiments 

2.  Experiments  on  the  strength  of  rocks  and  minerals 

3.  Static  pressure-volume  experiments 

4.  Dynamic  elastic  and  anelastic  experiments 

5 .  Thermodynamic  measurements 

6.  Consideration  of  the  form  of  equations  of  state 


The  contributions  from  these  sources  are  outlined  below: 


1 •  Shock  wave  experiments 

The  type  of  shock  wave  data  most  available  for  solids  including 
rocks  and  minerals  are  Hugonlot  curves.  The  Hugonlot*,  which  Is 
usually  depleted  In  the  pressure-volume  plane,  represents  the 
locus  of  states  that  a  material  may  achieve  upon  shock  compression 
from  some  Initial  state.  This  locus  represents  a  property  of  the 
SMterial  and  is  not  intrinsically  related  to  the  technique  employed 
in  achieving  the  shock  states.  The  Hugonlot  curve  which  is  most 
often  measured  is  termed  the  principal  Hugonlot;  this  corresponds 
to  an  initial  material  state  at  normal  temperature  and  pressure. 
However,  Hugonlot  curves  other  than  the  principal  Hugonlot  can  be, 
and  have  been, measured. 

Hugonlot  data  are  obtained  in  a  series  of  experiments  in 
which  plane  shock  waves  (for  which  the  pressure  profiles  are 
approximated  by  step  functions)  are  driven  into  a  series  of  material 
samples.  These  plane  shock  waves  produce  dynamic  one-dimensional 
compression  in  the  solid.  The  quantities  usually  measured  in 
such  experiments,  either  by  optical  or  electronic  techniques,  are 
the  shock  and  particle  velocities.  These  velocities,  and  the 
initial  density,  when  introduced  into  the  Rankine-Hugoniot  conservation 
equations  give  the  shock  pressure,  density,  and  change  in  internal 

_ 

Often  called  the  dynamic  adiabat  in  the  Soviet  literature 
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energy  produced  by  the  shock  wave.  This  type  of  experiment  is 
repeated  using  varying  incident  shock  pressures,  thereby  tracing 
out  the  entire  Hugonlot  curve.  The  Hugonlot  curve  then  represents 
a  pressure-density  curve  along  which  the  internal  energy  varies 
in  a  known  way.  This  is  somewhat  different  thar  the  case  of  an 
isotropic  or  an  Isothermal  pressure-density  curve  for  which  the 
entropy  and  temperature  are  fixed,  respectively. 

At  least  one  shock  pressure  range  exists  in  all  rocks  and 
minerals  in  which  the  final  high  pressure  shock  state  is  achieved 
by  two  (and  in  some  cases  more)  shock  fronts.  Multiple  shock 
fronts  arise  from  dynamic  yielding  or  shock-induced  polymorphism. 

In  the  case  of  dynamic  yielding,  the  amplitude  of  the  resulting 
shock  front  is  called  the  Hugonlot  elastic  limit  (HEL) .  The  HEL 
represents  the  maximum  normal  stress  that  the  material  can  withstand 
under  one-dimensional  compression  without  Internal  rearrangement 
taking  place  at  the  shock  front.  The  HEL  is  thus  an  important 
measure  of  the  stress  which  the  material  may  withstand  under  shock 
compression.  It  therefore  serves  as  a  useful  parameter  in 
formulation  of  an  equation  of  state  which  takes  into  account  the 
strength  of  the  material.  HEL  data  have  been  :ompared  to  static 
and  quasi-static  strength  data  in  order  to  formulate  a  complete 
yield  criteria. 

Hugonlot  data  for  rocks  and  minerals  are  often  tabulated  in 
terms  of  the  quantities  measured  in  the  experiments,  shock  and 
particle  velocity,  as  well  as  in  terms  of  pressure  and  density;  the 
latter  are  useful  in  formulating  complete  equations  of  state. 


Pressure,  kb 


Hugonlot  data  reported  to  1963  which  Included  a  few  rocks  are 
tabulated  by  Rinehart  (196$).  A  more  recent  and  complete  tabulation; 
which  Includes  recent  Soviet  and  Los  Alamos  results,  Is  given  by 
van  Thiel  (1966).  Data  gathered  by  McQueen  and  co-workers  for 
a  large  number  of  non-porous  rocks  and  minerals  within  the 
hydrodynamic  regime  is  collected  in  Birch  (1966)  and  some  are 
described  In  detail  by  McQueen  et  al  (1967b).  Additional  data 
In  this  regime  are  available  from  the  work  of  Trunin  e£  al_  (1965) , 
Al'tshulcr  (196$)  and  Jones  e£  al  (1968).  Hugonlot  data  for 
various  rocks  and  minerals  In  the  dynamic  yielding  regime  are 
reported  by  Ahrens  £t  £l  (1966),  Ahrens  (1969b);  and  Peterson 
it  il  (1969).  Hugonlot  data  for  porous  rocks  of  Interest  to  the 
Ferrls-Wheel  program  are  susecarlzed  In  Flanagan  (1966).  Shock 
measurements  on  various  porous  geologic  materials  performed  at 
Sandla  Corporation  and  Los  Alamos  have  been  summarized  by  Bass 
(1966)  and  McQueen  (1968).  Our  analysis  of  the  shock-wave 
data  has  largely  concentrated  on  developing  equations  of  state  at 
high-pressures  for  non-porous  rocks  and  minerals. 

Recently,  experimental  data  delineating  the  release  adiabats 
for  high-pressure  shock  states  of  several  geologic  materials  have 
become  available.  The  release  adlabat  curves  represent  the 
succession  of  states  through  which  a  mass  element  brought  to  high 
pressure  by  a  shock  front,  l.e.  shocked  to  a  Hugonlot  state, 
passes  through  upon  being  returned  to  ambient  pressure  and  achieving 
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a  post-shock  volume.  Hence,  a  family  of  release  adiabat  curves  exists 
for  each  Hugoniot  curve. 

A  limited  number  of  release  adiabat  data  have  been  reported  by 
Weidermann  (1964),  Anderson  £t  al  (1966),  for  porous  materials  and 
for  solid  materials  by  Ahrens  (1968),  Ahrens  e£  aJL  (1968), 

Rosenberg  £t  al  (1968),  and  Peterson  ££  til  (1969).  Although  this 
type  of  data  is  somewhat  more  difficult  to  obtain  experimentally, 
it  is  of  considerable  interest  from  materials  which  upon  shock 
compression,  either  undergo  irreversible  mechanical  crushing  or 
irreversible  phase  changes  accompanied  by  large  changes  in  volume. 

2.  Experiments  on  the  Strength  of  Rocks  and  Minerals 

A  considerable  body  of  data  relating  to  the  mechanical  strength 
properties  of  rocks  as  a  function  of  rock  structure,  stress  state, 
strain  rate,  and  moisture  content,  has  recently  been  summarized 
by  Walsh  and  Brace  (1969).  Earlier  tabulations  of  the  data  in  this 
area  are  available  through  the  AFCRL  Bibliography  of  Rock  Deformation 
(Riecker  et  al,  1965,  1966).  Data  on  the  strength  of  a  wide  class 
of  rocks  and  minerals  under  stresses  up  to  ~  40  kb  have  been 
summarized  by  Handin  (1966).  The  work  of  Brown  £t  jil  (1967) 
who  have  obtained  stress-strain  curves  from  various  rocks  under 
one-dimensional  static  compression  —  the  same  geometry  as  in 
shock  experiments  —  is  also  of  special  interest.  In  addition, 
we  have  found  the  recent  work  of  Towle  and  Riecker  (1968)  on 
the  shear  strength  of  minerals  at  high  temperature  and  pressure 
useful  in  reducing  the  shock-wave  data  (See  Section  III-2). 


3.  Static  Pressure  -  Volume  Measurements 


We  have  found  that  comparison  of  high-pressure  volume  data 
with  Hugoniot  data  is  useful  for  evaluating  the  effect  of  material 
strength  under  shock  loading  conditions.  Isothermal-pressure- 
volume  data  for  series  of  porous  and  solid  rocks  and  minerals 
have  been  obtained  with  piston-cylinder  type  apparatus.  These 
have  been  reported  by  Stephens  (1964),  Stephens  and  Lilley  (1966),  LaMori 
(1967)  and  La  Mori  (1968).  The  earlier  data  of  Bridgman  (tabulated  in 
Birch,  1966)  most  of  which  is  below  30  kb,  is  in  several  Instances 
the  only  static  data  available  for  a  given  mineral.  Drickamer 
(Perez-Albuerne ,  and  Drickamer,  1965,  Lewis  and  Drickamer,  1966) 
and  McWhan  (1967)  have  obtained  an  Important  series  of  pressure- 
volume  curves  for  oxides  to  approximately  40  kb  using  high 
pressure  X-ray  techniques.  More  recently  the  work  of  the  high- 
pressure  X-ray  group  at  Rochester  (Bassett  &£  al,  1968;  Bassett 
and  Barnett,  1959;  Mao  et  al,  1969;  Liu  e£  al ,  1969)  have 
significantly  Increased  our  knowledge  of  the  pressure-volume  behavior 
of  high  pressure  polymorphs  of  the  silicates. 

4.  Dynamic  Elastic  and  Anelastic  Experiments 

The  measurements  of  elastic  wave  velocities  and  elastic 
moduli  of  rocks  and  minerals  are  summarized  by  Birch  (1966),  and 
Anderson  and  Liebermann  (1966) ,  and  for  a  wide  range  of  compounds 
by  Simmons  (1965)  and  for  oxides,  Anderson  et  al  (1968).  Measurements 
of  elastic  moduli  to  pressures  of  10  kb  at  temperatures  of  800°C 
are  now  possible  and  have  been  carried  out  for  selective  materials 
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which  may  be  used  with  ultrasonic  measuring  techniques  (see  for 
example,  Spetzler,  1969  a,  b). 

The  attenuation  of  stress  waves  In  the  seismic  or  so-called 
elastic  regime  Is  controlled  by  the  small  intrinsic  anelastlc 
properties  of  geologic  materials  as  well  as  the  extrinsic  loss 
of  energy  due  to  Internal  interfaces  such  as  joints,  faults, 
and  bedding  surfaces  which  occur  in  virtually  all  crustal  rocks. 

Data  for  various  engineering  materials  and  for  rocks  obtained 
from  laboratory  measurements  to  about  1964  are  reviewed  by  Knopoff 
(1965a).  Additional  data  particularly  for  porous  and  sedimentary 
rocks  are  Included  In  a  data  compilation  given  by  Bradley  and 
Fort  (1966).  Various  anelastlc  mechanisms  in  the  earth  have  been 
discussed  by  Savage  (1964),  Gordon  and  Nelson  (1966),  and 
Merkulova  (1966).  More  recently,  a  review  of  pertinent  anelastlc 
data  for  ceramics  and  geologic  materials  has  been  prepared  by 
Jackson  and  Anderson  (1969), 

5,  Thermodynamic  Measurements 

The  specification  of  states  which  are  displaced  from  the 
normally  measured  pressure-volume  curves  In  the  pressure-volume 
plane,  such  as  the  principal  Hugoniot,  requires  supplementary 
data  such  as  the  Grunelsen  ratios  and  the  volume  expansion  coefficient. 
Also,  shock  temperature  calculations  require  knowledge  of  specific 
heat,  usually  assumed  to  be  a  function  of  volume  and  temperature. 

In  order  to  extend  the  equation  of  state  through  a  regime  where  a 
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solid-solid  transition,  melting  or  perhaps  vaporization  takes  place, 
supplementary  data  such  as  the  enthalpy  of  the  appropriate  species 
ece  required. 

Thermodynamic  data  for  minerals  at  standard  conditions  are 
given  by  Roble  (1966  and  1968).  Thermal  expansion  data  for 
minerals  are  tabulated  by  Skinner  (1966) .  High  temperature  enthalpies 
for  many  of  the  species  found  In  earth  materials  are  available  In 
the  Janaf  Tables  (Stull  £t  al,  1965).  Debye-temperatures  and  the 
acoustic  Grunelsen's  ratio  are  extrapolated  from  ultrasonic  data 
such  as  summarized  by  Anderson  and  Liebermann  (1966),  Anderson 
(1968),  and  Anderson  et  al  (1969). 

6.  Equations  of  State 

Hugonlot  and  other  forms  of  equation  of  state  are  usually 
formulated  In  terms  of  expressions  containing  1,  2,  or  more 
adjustable  parameters.  These  relate  pressure  to  density  along  a 
particular  curve  (Hugonlot)  or  a  thermodynamic  path  (adlabat  or 
Isotherm).  States  along  different  pressure-volume  curves,  such 
as,  for  example,  two  states  differing  in  internal  energy,  E,  and 
pressure,  P,  at  the  same,  V,  are  often  related  by  the  Mie-Grunelsen 
equation 

P  -  PK  -  ^  (E-E^  (1) 
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P  and  Ev  are  the  pressure  and  energy  along  some  reference  curve. 
Here  y  is  the  Grunelsen  ratio  which  may,  at  high  temperatures, 
be  Independent  of  temperature.  (For  silicates, y  •  (3P/9E)y,  at 
around  room  temperature,  probably  depends  on  temperature  as  well 
as  volume.)  For  metals  it  is  often  assumed  that  y  depends  only 
on  volume.  Among  the  more  often  used  equations  of  state  are  the 
Mumaghan  equation 

P  -  <K  /n)[x3n-l]  (2) 

o 

and  the  Blrch-Mumaghan  second  order  equation 

P  -  A(X7-X5)  -  B  (X7-X5) (X2-l)  (3) 

which  as  Thomson  (1969)  has  pointed  out,  corresponds  to 
assuming  a  Eulerian  form  for  the  strain  in  finite  strain  theory. 
Here  Kq  and  n,  and  A  and  B  are  adjustable  parameters  and 

x  -  (p/p  )1/3  “  (V  /V)1/3  (4) 

o  o 

where  pq  and  p  are  the  initial  and  final  density  and  Vq  and  V 
are  the  Initial  and  final  volume.  Equation  3  may  be  written 


as 
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3K 

P - 2^  [X7-X5] tl-C(X2-l) ]  (5) 

where  K  ■  2A/3  is  the  bulk  modulus,  5  ■  B/A  and  K*  ■  (dK  /dP)  ■ 
o  o  o  o 

4 

4  -  y  5.  Thomsen  has  recently  discussed  an  alternate  form  of 
equation  5  corresponding  to  expressing  the  strain  In  the  theory 
of  finite  strain  In  Lagrangian  form.  The  equation  which  results 
Is 


P  “  I  Ko  t™"1]  f  Ko 


(6) 


Equation  6  gives  a  similar  pres sure- volume  curve  as  Equation  5 
for  compressions,  p/p  of  less  than  1.2.  Some  difficulties  are, 

O) 

however,  noted  In  applying  this  equation  to  pressure  volume 
data;  these  are  described  In  Section  III-4.  In  the  above 
equations  the  zero-pressure  bulk  modulus,  Kq  is  related  to  a  bulk 
sound  speed }  Cq jSnd  seismic  parameter,  $ ,  by 


K  /p  -  C  -  $ 
o  o  o 


(7) 


As  Indicated,  Is  the  pressure  derivative  of  the  bulk  modulus. 

In  practice  It  Is  simpler  to  use  the  Blrch-Mumaghan  equation  to 
fit  a  series  of  P-V  data  or  a  series  of  calculated  P-V  points 
than  it  Is  to  use  the  Murnaghan  equation  because  In  equation  2 
the  exponents  of  X  are  fixed  and  hence  evaluation  of  the  parameters 
A  and  B  Involves  only  the  solution  of  a  series  of  linear  simultaneous 
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equations . 

It  should  be  noted  that  both  the  Birch-Murnaghan  and  Murnaghan 

equations  place  restraints  on  the  high-order  parameter  9K/3P 

along  the  curve  which  they  define  in  the  P-V  plane.  The 

Murnaghan  equation  9K/9P  is  a  constant  while  in  the  Birch-Murnaghan 

equation  the  derivative  92K/9P2  is  fixed  by  the  form  of  equation  5. 

Another  two-or  three- parameter  equation  may  be  obtained  from 

the  often-used  linear  or  quadratic  shock  (U  )  and  particle  velocity 

s 

(U  )  relation 
P 

V  -  C  +  S]U  +  S2U2  (8) 

so  p  p 

where  Cq  represents  the  bulk  sound  speed  and  Si  and  S2  are  equation 
of  state  parameters.  In  a  case  of  a  two-parameter,  linear  equation 
(S2  ■  0)  it  follows  from  the  Rankine-Hugonlot  shock  equations  that 

P  -  C2  (V  -V) / IV  -S j  (V  -V) 2 ]  (9) 

OO  0  0 

Equation  9  is  particularly  well-suited  for  fitting  raw  Hugoniot 
data.  However,  the  thermodynamic  strains  posed  by  the  linear  U 

s 

and  Up  relation,  which  have  been  recently  discussed  by  Ruoff 
(1967)  reduces  the  usefulness  of  this  relation  by  forms  of  equations 
of  state  other  than  the  Hugoniot.  For  the  three-parameter  Aquation 


where  Si  is  either  obtained  experimentally  from  shock-wave  experiments 
or  from  ultrasonic  measurements  using  the  relation 
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(14) 


Recently,  Thomsen  (196.)  has  foliated  a  three-parameter  a, nation  of 
atata  alao  basal  on  flnlta  attain  thaory  which  traata  tha  Gronaiaan 

parameter  Y»  («/»«„.  and  3W  l„  .  aalf-conalatant  way.  Thla 
three-parameter  relation  is  given  by 

p  -  3K  a/a)  [n  -f  n2r  +f  „3AJ  +  X 

r  (15) 

|7+3n[l-^  (1-  T  cv/et>]  | 

where 


n  -  \  [X-2-1J 


a^ 


ET(/Cv) 


A  -  K  K'  ’  -  Ck1  -  -)  r  -  A  iL 
0  1  o  3'  r  9  7~ 


- 


yV-  2  a 


+  2v^ 


1  +  T 


cvc/eto 


)] 


(16) 


(17) 
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(  3T  )  pa0  ETO/CVO 


(18) 


Y  * 


r  - 


(19) 


(20) 


Here  a  Is  a  measure  of  a  length  the  unit  cell,  a  Is  the  volume 
expansion  coefficient,  K  Is  the  bulk  modulus  and  the  primes 
Indicate  differentiation  with  respect  to  pressure.  Also  E  is 
the  thermal  vibration  energy  and  Cy  Is  the  specific  heat  at  constant 
volume, the  subscripts  o  indicate  standard  conditions.  Because 
the  Mie-Gruneisen  assumption,  y  ■  y(V),  is  implicit  in  equations 
14  through  20,  equation  15  is  expected  to  be  useful  for  geologic 
materials  at  or  above  Debye-temperature  of  the  material.  This 
equation,  therefore,  has  considerable  promise  in  describing  these 
materials  at  high  temperatures  and  in  the  hydrodynamic  regime. 

III.  Analysis  of  Equation  of  State  Data 

The  major  emphasis  in  research  carried  out  during  the  past  year 
has  been  the  analysis  of  a  large  body  of  available  shock-wave  data 
for  rocks,  minerals,  and  earth  forming  oxides.  In  the  analysis  of  the 
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data  we  have  attempted  to  use  the  physically  significant  parameters 
and  incorporate  supplementary  data  such  as  crystal  structure,  atomic 
weight,  elastic  moduli  at  normal  and  elevated  temperatures  and 
pressures,  as  well  as  thermodynamic  properties  and  strength  of  material 
data.  The  research  performed  on  this  contract  has  been  presented 
in  detail  in  five  manuscripts  which  are  in  various  stages  of  publication 
as  indicated  below. 

1.  Equation  of  state  and  crystal  structure  of  high-pressure 

phases  of  shock  silicates  and  oxides,  Thomas  J.  Ahrens,  Don  L.  Anderson, 
and  A.  E.  Ringwood,  Reviews  of  Geophysics ,  7_,  667-708,  1969. 

2.  Equation  of  state  and  pressure  dependence  of  the  shear  strength 
of  two  oxides  under  shock  compression,  Geoffrey  Davies  and 
Thomas  J.  Ahrens,  to  be  published. 

3.  The  Equation  of  state  of  stishovite,  Thomas  J.  Ahrens, 

T.  Takahashi,  and  Geoffrey  Davies,  J.  Geophvs.  Res.,  in  press,  1970. 

4.  Calculation  of  Gruneisen's  ratio  of  Periclase  to  10  kb  and 
1200°K  from  Ultrasonic  Data,  Thomas  J.  Ahrens,  and  Hartmut  A.  W. 
Spetzler,  to  be  published. 

2- 

5.  Heat  of  Formation  of  0  ,  Edward  S.  Gaffney  and  Thomas  J.  Ahrens, 

J.  Chem,  Phys.,  51,  1088-1091,  1969. 

The  results  obtained  in  these  papers  are  outlined  and  summarized 


below: 
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1.  Equation  of  State  and  Crystal  Structures  of  High  Pressure 
Phases  of  Shocked  Silicates  and  Oxides 

The  raw  tlugonlot  data  for  minerals  (Figure  1)  may  often  be 
represented  as  lylnR  In  one  or  more  of  the  three  regimes  which  are 
designated  as  the  low  pressure  phase,  mixed  phase,  and  high  pressure 
phase  (hpp)  regime. 

In  this  paper  the  shock  wave  data  for  rocks  and  minerals 
obtained  by  McQueen  and  Marsh  (1966),  Me  Queen  e_t  al  (1967b), 

Al' shuler  c^t  a_l  (1965),  Trunin  c_t  aL  (1965) ,  Wackerle  (1962), 
and  Ahrens  c_t  al  (1967b)  were  analyzed  In  order  to  obtain  equation- 
cf-statc  properties  and  Infer,  especially  from  the  zero-pressure 
density,  the  crystallography  of  the  high  pressure  phase  (hpp) . 

With  the  exception  of  stlshovlte  (S10,)  there  Is  no  direct 
knowledge  of  the  crystallography  of  the  hpp  produced  In  the 
shock  experiments. 

For  manv  of  the  subject  materials  the  low-pressure  regime 
corresponds  to  stress-voluaie-encrgy  states  In  material  that  has 
yielded  under  one-dlmenslonal  compression  and  has  thus  been  shocked 
Into  states  which  are  In  an  effectively  hydrostatic  state  of  stress, 
l.e.,  the  principal  stresses  arc  nearly  equal.  In  some  cases, 
however,  for  exasple,  AljOj  and  MgO,  significant  stress  differences 
appear  to  be  present  at  Hugonlot  states  up  to  substantial  stress 
levels.  These  are  discussed  below.  The  mixed  phase  region  is 
assumed  to  represent  a  mixture  of  both  the  high  and  low  pressure 
phase  material  while  states  In  the  high-pressure  regime  are 
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assumed  to  represent  material  that  has  been  wholly  converted  to 
the  high-pressure  phase.  In  this  paper  we 

analyzed  the  shock  data  in  the  high  pressure  region.  Although 
in  practice  it  is  often  difficult  to  ascertain  the  lower  pressure 
limit  of  this  regime.  It  should  be  noted  that  the  states  lying 
in  the  high  pressure  regime  are  achieved  by  direct  shock  transition 
from  the  low  pressure,  initial  phase  which  is  at  standard  conditions. 
Thus,  no  assumption  of  the  equation  of  state  of  the  low  pressure 
phase  material  need  be  made. 

The  parameters  in  the  Birch-Mumaghan  equation  were  obtained 
for  the  24  materials  analyzed  along  the  raw  Hugoniot  centered  at 
zero  pressure,  the  metastable  Hugoniot  of  the  shock-induced  hpp 
and  the  adiabat  and  25°C-isotherm  of  the  hpp  (Table  1  and  2). 

Of  all  the  minerals  and  rocks  for  which  data  are  analyzed,  alx 
but  three  (A^O^  -  corundum,  MgO  -  periclase,  and  Mn02  -  pyrolusite) 
appear  to  have  at  least  one  major  phase  transition  between  0  and 
800  kb.  Because  of  the  marked  decreases  in  volume  (typically 
20-50%),  the  zero-pressure  density  of  the  hpp*s  were  obtained 
by  restraining  the  hpp  adiabat,  calculated  from  the  Hugoniot  data, 
so  that  the  zero-pressure  density,  pq,  zero-pressure  seismic 
parameter,  pq,  and  mean  atomic  weight,  M,  satisfied  the  equation 

p  /M  =>  A  $n 
o  o 


(20) 
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suggested  by  Anderson  (1967),  here  A  *■  0.048  and  n  ■  0.323. 

(A  recent  revision  of  his  analysis  (Anderson,  1969)  suggests 

values  of  A  ■  0.0492  and  n  ■  1/3  are  more  appropriate.  The  new 

value  of  A  and  n  give  values  of  zero-pressure  density  which  are 

1  to  4%  lower  than  those  reported  here.)  When  this  restraint  is 

applied  to  the  shock  data  for  corundum,  periclase,  pyrolusite, 

and  stishovite,  for  which  the  zero-pressure  derivatives  are 

known,  the  error  in  the  inferred  zero-pressure  density  varies 

from  -0.82  for  stishovite  to  +2.82  for  MgO.  The  transformation 

energy  Efcr,  at  standard  pressure  and  temperature  is  also  employed 

in  the  analysis.  In  general,  transformation  energies  are  estimated 

from  Independent  thermochemical  data  for  the  same  or  analogous 

transitions.  Uncertainties  of  •-  50%  in  E_  result  in  uncertainties 

tr 

of  only  ~  12  in  the  calculated  zero-pressure  densities.  The 
sensitivity  of  calculated  equation-of-state  parameters  to  assumptions 
which  are  made  concerning  both  the  initial  (zero-pressure)  value 
of  Gruneisen's  ratio  (y) ,  and  the  dependence  of  y  on  compression 
were  tested.  Since  the  initial  value  of  y  can  vary  by  almost  a 
factor  of  2,  depending  upon  whether  the  Slater  or  the  Dugdale- 
McDonald  value  is  used,  it  is  found  that  the  assumption  of  the 
initial  value  of  y  affects  the  calculated  equation-of-state  to 
a  much  greater  degree  than  the  assumed  dependence  of  y  on 
compression.  For  example,  in  Fe^O^  a  change  in  y  fr^m  0.34  to 
0.63  (Slater  to  Dugdale-McDonald)  changes  the  inferred  zero-pressure 
density  of  the  hpp  by  only  ~  0.32  and  the  zero-pressure  bulk  modulus 
by  ~  1.5%. 
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The  probable  crystallography  of  the  oxide  and  silicate 
phases  produced  at  high  shock  pressures ,  is  inferred  from  the 
calculated  zero-pressure  density  and  indirect  lines  of  evidence, 
such  as  high-pressure  petrologic  studies  of  analog  systems 
(e.g.,  germanates,  titanates)  under  static  conditions  as  well  as 
the  classical  rules  of  crystal  chemistry*  Except  for 
stishovite,  shock  recovery  experiments  have  not  as  yet  yielded 
samples  of  the  hpp * s  that  are  actually  produced  in  shock  experiments. 
Statements  as  to  the  probable  crystallography  of  the  hpp's  are 
therefore  necessarily  subject  to  considerable  uncertainty.  This 
uncertainty  largely  depends  on  the  availability  of  data  for 
appropriate  analog  systems.  The  oxides,  A^O^,  MgO,  and  Mn02 ,  do 
not  appear  to  have  significant  phase  transitions  up  to  ~  1200  kb. 

Rutile,  Ti02>  appears  to  transform  from  an  initial  density  of 

3  3 

4.25  g/cm  to  a  hpp  with  a  zero-pressure  density  of  5.71  g/cm  . 

As  previously  inferred  by  McQueen  ££  al. (1967a)  this  hpp  probably 

has  the  fluorite  structure.  Although  similar  transformations 

are  known  to  occur  in  Ge02  and  Pb02,  the  analogous  transformations 

are  not  detected  in  the  analysis  of  the  Si02  shock  data  (available 

to  2000  kb  and  presumably  representing  stishovite  above  ~  400  kb). 

3 

The  hpp  of  hematite,  Fe^^,  density  5,96,  vs.  5.27  g/cm  for  the 
low  pressure  phase,  is  inferred  to  have  the  perovskite  structure 
although  the  B-rare  earth  structure  (SCjO^)  is  another  possibility. 

Spinel,  MgAljO^,  transforms  to  a  phase  which  has  zero-pressure 

3  3 

density  of  4.19  g/cm  ,  which  compares  with  a  density  of  3.86  g/cm 

for  the  isochemical  mixed  oxides.  We  infer  a  structure  for  this 
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phase  similar  to  that  of  CaFe^^.  The  hpp  of  magnetite,  Fe.^,  has 

3 

a  density  of  5.54  as  compared  with  5.20  g/cm  for  the  low-pressure 

phase.  This  phase  may  represent  either  a  disproportionated 

product  of  FeO  +  Fe^^  (perovskite)  or  possibly  the  orthorhombic 

CaFe.O,  structure. 

2  4 

The  minerals,  albite,  oligoclase,  and  microcline  all  appear 

3 

to  transform  to  the  hollandlte  structure,  density  ~  3.85  g/cm 
(Ringwood  et  al.  (1969).  Although  this  transformation  does  not 
appear  to  be  completely  over  the  range  of  data  in  the  pure  mineral 

samples,  the  zero-pressure  density  of  the  high-pressure  products 

3  3 

of  Westerly  granite  (3.96  g/cm  )  and  anorthosite  (3.51  g/cm  ) 

are  consistent  with  a  mixture  of  stishovite  and  feldspar  (in 

the  hollandlte  structure)  existing  at  high  pressure. 

In  the  case  of  olivines,  since  the  hpp  zero-pressure  density 

is  typically  some  20%  higher  than  the  starting  material,  the 

olivine-spinel  transition  (usually  involving  a  density  increase 

of  less  than  10Z)  is  clearly  overdriven.  Although  the  hpp  densities 

are  close  to  those  which  would  occur  if  this  mineral  were  to 

disproportionate  either  to  the  oxides-periclase ,  wustite  and 

stishovite-or  to  ilmenite  and  rock  salt  structure,  in  the  range 

Fogg  to  Fo^g,  it  appears  most  likely  that  the  SrjPbO^  structure 

[Wadsley  et  al,  1968]  is  produced.  In  the  case  of  pure  forsterite, 

3 

with  an  initial  density  of  3.21  g/cm  ,  a  hpp  material  having  a 

3  3 

density  of  4.31  g/cm  is  obtained.  This  compares  to  4.12  g/cm 

for  the  hpp  of  Twin  Sisters  dunite  (Fogg).  A  somewhat  higher 

density  of  the  magnesia  end -member  may  represent  the  formation 

of  the  K2NiF^  rather  than  the  Sr2PbO^  structure.  The 


K2NiF4 


structure  is  taken  by  the  olivines  Ca^GeO^  and  MpCaGeO^  (Ringwood 
and  Reid,  1969).  This  result  should,  however,  be  taken  with 
caution  as  only  a  few  less  reliable  Hugoniot  points  are  reported 
for  this  material. 

The  Hugoniot  data  for  bronzitite  (En^g)  gives  a  hpp  zero- 

3 

pressure  density  of  3.74  g/cm  as  compared  to  a  starting  density 
3 

of  3.27  g/cm  .  The  hpp  is  thought  to  represent  a  garnet  or 
ilmenite  structure  or  possibly  a  disproportionation  into  spinel 
plus  stishovite.  The  relatively  low  zero-pressure  bulk  modulus 
inferred  for  the  hpp  material  (~  2100  kb)  suggests  that  the 
experimental  data  may  not  all  lie  within  the  high  pressure  phase 
regime. 

As  in  the  case  of  pure  forsterite,  the  data  for  pure  enstatite 

are  also  less  reliable.  These  also  give  significantly  higher 

densities  than  the  material  containing  more  FeO.  The  inferred 

3 

zero-pressure  density  of  the  hpp,  4.20  g/cm  ,  correlates  well 

3 

with  the  expected  zero-pressure  density,  of  4.25  g/cm  for  the 
perovskite  structure. 

Sillimanite  and  andalusite  have  hpp  zero-pressure  densities 
3 

of  4.00  and  3.95  g/cm  .  These  values  probably  represent  the 
disproportionation  reaction  products,  A^O^  and  Si02  (stishovite). 
This  mixture  would  have  a  theoretical  zero-pressure  density  of 
4.09  g/cm"*. 

The  equation-of-state  parameters  of  the  hpp's  display  several 
trends  which  are  believed  significant.  We  find  that  the 


second  order  parameter  £  in  the  Birch-Murnaghan  equation  is 

nearly  zero  for  the  oxides  MgO  and  AljO^.  For  stishovite,  we 

find  £  ■  0.73  and  K'  ■  (dK/dP)  *  3  (a  later  analysis  yields 

os  s 

a  value  of  K'  of  7,  Ahrens  et  al,  1969a)  while  for  the  hpp  of 
os  —  — 

olivines  and  pyroxenes  values  of  £  closer  to  1  with  values  of 

K'  of  2  to  3  are  obtained.  The  value  of  dK/dP  calculated  along 
os 

the  adiabat  and  Hugoniot  curves  (the  latter  representing  higher 
temperatures) ,  decrease  with  increasing  compression.  The  rate 
of  decrease  of  K'  with  density  inferred  from  the  Birch-Murnaghan 
equation  is  similar  for  all  the  oxides  and  hpp's  of  the  silicates. 

3 

These  vary  from  -0.5  to  -0.7  cm  /g.  Increasing  the  density  upon 

substitution  of  FeO  for  MgO  gives  values  for  the  change  of  (dK/dP) 

3  3 

with  density  in  the  range  -0.3  cm  /g  to  -1.6  cm  /g.  The  values 
of  y  for  the  hpp's  are  between  0.5  and  1.0.  These  are 
related  to  which  is  also  generally  lower  for  the  hpp's  and 
decrease  with  compression  and  increasing  iron  content. 

Since  the  Hugoniot  data  for  MgO  and  AljO^  indicate  that  no 
significant  phase  transition  has  taken  place  in  these  materials, 
the  zero-pressure  densities  were  constrained  to  their  zero-pressure 
values  for  the  purpose  of  further  analysis.  The  resulting  bulk 
moduli  obtained  were  higher  than  the  ultrasonic  values  of  1622 
and  2521  kb.  Since  no  internal  rearrangement  takes  place  at 
the  shock  front  (presumably  due  to  the  lack  of  a  phase  transition) , 
strength  effects  will  produce  a  Hugoniot  curve  for  MgO,  AljO^,  and 
MnOj,  which  is  offset  in  shock  stress  above  the  hydrostatic  curve. 
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Upon  correcting  the  raw  Hugoniot  data  for  MgO,  ceramic  and  crystalline 
A^O^,  for  residual  strength  effects  (by  relating  the  net  shear 
stress  at  high  shock  pressures,  to  that  existing  at  the  Hugoniot 
elastic  limit  (Ahrens  e£  £l  ,  1968a))  an  analysis  of  the  hydrostatic 
Hugoniot  was  performed.  A  corrected  Hugoniot  yields  values  for  the 
zero-pressure  adiabatic  bulk  modulus  of  1675,  2563,  and  2919  kb  for 
MgO,  ceramic, and  crystalline  A^O^.  These  are  in  good  agreement 
with  the  ultrasonic  data.  (A  more  recent  analysis  of  the  shock 
data  in  which  the  shear  strength  is  allowed  to  increase  with  mean 
principal  stress  gives  similar  results  (Davies  and  Ahrens,  1969).) 

The  good  agreement  in  both  bulk  moduli  and  (dKQs/dP)s  is  consistent 
with  the  result  that  MgO  and  A^O^  not  unt^er80  phase  changes 
over  the  range  of  the  shock  data  and  furthermore  that  these 
materials  support  significant  levels  of  shear  stress  above  shock 
pressures  of  150  and  300  kb,  respectively. 

2.  Equation  of  State  and  Pressure  Dependence  of  Shear  Strength 
of  Two  Oxides  under  Shock  Compression 

Isotherms  and  a  family  of  release  adiabats  are  calculated 
from  the  periclase  (MgO)  and  corundum  (A^O^)  Hugoniot  data.  A 
strength  correction  which  depends  on  pressure  and  temperature 
in  accordance  with  the  empirical  relation  proposed  by  Towle  and 
Riecker  (1968),  is  used  to  bring  the  zero-pressure  and  high- 
temperature  densities  into  self-consistent  agreement  with  the 
high-temperature  thermal  expansion  and  ultrasonic  data  (Table  3.). 


Assuming  a  linear  dependence  of  shear  strength,  S,  with  pressure 
according  to 


S  -  S  •  (1  +  cP)  (21) 

o 

as  suggested  by  Towle  and  Riecker,  has  improved  agreement  of  the 
shock  data  for  single  crystal  corundum  with  the  zero-pressure 
and  high  temperature  ultrasonic  data.  Here  S  is  the  strength  at 
hydrostatic  pressure  P,  Sq  is  the  initial  strength  and  e  is  a 
constant.  For  a  periclase,  the  agreement  between  the  high  temperature 
parameters  of  the  adiabats  calculated  and  those  observed  has 
also  been  improved  by  the  amount  comparable  to  that  obtained  earlier 
(Ahrens  ejt  al,  1959c)  by  assuming  a  constant  strength  which  is 
independent  of  pressure.  In  this  case,  errors  in  the  analysis  and 
in  the  data  may  be  of  the  same  order  as  the  effects  obtained.  For 
ceramic  corundum  this  agreement  increases  with  increasing 
temperature  and  cannot  be  substantially  reduced  using  this  method. 

The  pressure  dependence  of  a  shear  strength  obtained  from  this 
analysis  of  shock-wave  data  is  indicated  in  Table  4  and  is  in 
surprisingly  good  agreement  with  the  pressure  dependence  of  the 
shear  modulus  which  has  been  hypothesized  by  Towle  and  Riecker.  In 
the  analysis  of  the  shock  wave  data,  we  have  again  assumed  the 
Birch-Murnaghan  form  (equation  5)  for  the  pressure-volume  curve. 

The  errors  introduced  by  employing  this  assumption  rather  than 
using  some  higher  order  equation  or  another  two-parameter  formula 
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are  difficult  to  determine.  One  possible  way  to  estimate  the 
errors  involved  is  to  compare  the  results  obtained  with  those 
obtained  by  other  equations  of  state  such  as  the  Murnaghan 
equation  or  equations  derived  from  different  interionic  potentials. 
Thomsen  (1969)  quoted  an  equation  derived  from  finite  strain  theory 
using  the  Lagrangian  definition  of  strain  rather  than  the  Eulerian 
definition.  This  two-parameter  equation  as  given  by  equation  6 
was  fitted  to  the  Hugoniot  data;  but  resulting  values  of  Kq  are 
in  poor  agreement  with  the  ultrasonic  data  and  the  values  of 
were  unrealistically  large,  i.e.,  10  to  15.  Thus,  this  equation 
appears  to  be  a  very  poor  representation  of  the  binding  forces 
involved.  The  accuracy  of  the  Birch-Murnaghan  equation  is  still 
not  determined,  but  because  of  its  success  representing  much  of 
the  shock-wave  data  probably  gives  Kq  to  within  about  5%.  As 
indicated  in  Table  4  the  pressure  dependence  of  the  shear  strength 
obtained  for  periclase  and  corundum  is  comparable  to  the  first- 
order  pressure  dependence  of  the  shear  modulus  and  supports  the 
hypothesis  of  Towle  and  Riecker  (1968). 

3.  Equation  of  State  of  Stishovite 

With  the  recent  availability  of  new  Hugoniot  data  for  sandstone 
(Jones  et  al,  1968),  fused  quartz  (McQueen,  1968;  Shipman  and  Isbell, 
1969)  we  have  reanalyzed  the  high-pressure  data  for  SiOj  in  the  high 
regime  where  quartz  has  transformed  to  stishovite.  Combining 
these  new  data  with  those  of  Wackerle  (1962)  and  Al'tshuler  (1965) 
for  quartz  (Figure  2),  and  using  the  procedure  outlined  by 


McQueen  ejt  al.  (1963)  an  Independent  determination  of  the  Gruneisen 
ratio  at  high  pressures  is  obtained  (Figure  3).  Using  this 
determination  of  y  and  the  recently  available  data  for  the  zero-pressure 
density,  transformation  energy  of  quartz  to  stishovite,  specific 
heat  and  thermal  expansion  coefficient  for  stishovite  an  isothermal 
and  adiabatic  equation  of  state  is  constructed  (Table  5). 
Stishovite  has  a  surprisingly  high  zero-pressure  Gruneisen  ratio 
of  approximately  1.7.  This  compares  to  ~  0.9  obtained  earlier 
by  McQueen  eit  al  (1963)  and  ourselves  (Ahrens  et  al,  1969c). 

We  find  that  the  value  0.9  is  consistent  with  the  earlier  calculations, 

but  only  at  higher  pressures.  It  appears  that  the  rapid  decrease  of  y 
along  the  Hugoniot  from  a  value  of  ~  1,7  at  standard  temperature 
and  pressure  to  0.9  kbj may  largely  be  a  temperature  effect. 

4.  Calculation  of  Gruneisen 's  ratio  of  Periclase  to  70  kb  and 
1200°K  from  Ultrasonic  data 

The  recent  availability  of  high-temperature  and  high-pressure 
ultrasonic  data  from  minerals  such  as  MgO,  has  permitted  a  cricital 
evaluation  of  thermodynamic  Gruneisen' s  ratio  as  a  function  of 
temperature  and  pressure.  In  this  analysis  we  apply  a  form  of  the 
Gruneisen  ratio 


y  -  V(3P/9E)V 


(22) 
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directly  to  the  ultrasonic,  thermal  expansion,  and  enthalpy  data. 

The  zero-pressure  adiabatic  bulk  modulus,  Kqs  and  its  pressure 

derivatives  (3K  / 3P)_  have  been  determined  ultrasonically  by 

os  I 

Spetzler  (1969b).  After  correcting  Spetzler's  mixed  derivatives 
to  adiabatic  conditions,  and  applying  the  thermal  expansion  data 
for  MgO,  a  series  of  pressure-volume  adiabats  were  constructed  using 
the  data  of  Victor  and  Douglas  (1958)  for  the  high-temperature 
enthalpy  of  MgO.  Gruneisen’s  ratio  was  then  calculated  as  a  function 
of  pressure  and  volume  (Figure  4). 

The  most  interesting  result  of  these  calculations  is  the 
marked  decrease  of  y  with  increasing  temperature  at  constant 
volume  (10%  at  normal  volume). 

The  often  employed  assumption  of  the  Mie-Gruneisen  equation  of 
state  which  implies  that  y  is  only  a  function  of  volume,  may  be 
in  error  below  the  Debye  temperature  (940°K) .  If  the  Mie-Gruneisen 
assumption  were  true  the  lines  of  constant  y  shown  in  Figure  2 
should  be  vertical.  As  can  be  seen  they  are  actually  nearly 
parallel  to  the  adiabats  (or  isotherms)  in  this  pressure  and 
temperature  range.  Because  of  the  increased  excitation  of  optical 
phonon  modes  with  temperature  between  300  and  1200°K,  the 
Gruneisen  ratio  of  MgO  appears  to  be  more  dependent  on  temperature 
than  volume  in  this  regime. 

According  to  the  phonon  dispersion  curves  given  by  Peckham 

12 

(1967)  the  maximum  acoustic  vibration  frequency  is  ~  15.5  x  10  cps. 


The  relation  T  •  h\*/k  give*  a  temperature  for  end  tat  Ion  of 
all  the  acoustic  node*  of  7))*.  The  average  acoustic  Crunelsen 
ratio  give  theoretical  value*  of  1.19  to  l.4|. 

The  resulting  Crunrl«en  ratio  for  the  transverse  and  longitudinal 
optic  node*  Rive  an  average  of  I . H  i  0.0 h.  If  at  744*K  only 
the  acoustic  node*  are  excited,  we  expect  >  •  >  •  1.40.  At 

In  AC 

high  temperature*  (above  I0W),  both  the  acoustic  and  optic 

Mode*  should  show  a  lover  composite  Crunelsen  ratio  given  by  i ^  *  1/2 

(  »  .  ♦  >  )  *  *»  l.y>.  The  value  of  the  thermodynamic  Crunelsen 

AC  Op 

ratio  which  la  measured  (~  I.W)  la  srtervhat  higher  than  the 
theoretical  estimate  of  1.40.  The  theoretical  fitlMte  doea 
however  Indicate  a  decrease  In  gaiwu  with  Inrreaalng  temperature. 

Ve  auggeat  thla  result  arise*  Iran  the  successive  excitation  of  the 
low  >  optical  mode*  with  Increasing  temperature.  Above  •»  1016* 
the  Crunelaen  ratio  should  become  nearly  Independent  of  temperature 
and  the  Kle-Crunel*en  equation  should  he  valid. 

J.  Heat  of  formation  of  0* 

An  part  of  our  continued  Interest  In  the  nature  of  the  high 
preaaure  phase*  produced  In  shocked  silicates  and  oxide*  a  series 
of  Born-Harer- type  calculation*  were  used  to  calculate  the  lattice 
energies  of  simple  oxides  fHgO,  BeO,  CaO,  and  7.nO) .  Repulsion 
and  other  non-Coulombl c  contribution*  to  the  lattice  energy  are 
obtained  using  thermodynamic  and  recent  ultrasonic  data  for  the 
bulk  moduli  and  the  Isothermal  pressure-temperature  derivatives 
of  the  elastic  constants.  Using  thermochemical  data  for  the 
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heat  of  formation  of  MgO,  CaO,  and  BeO  and  their  cations,  the 

2-  °  2- 

heat  of  formation  of  0  ,  AH^(0  ),  is  calculated  to  be  197  ±  5  kcal/ 

2- 

molc.  Using  the  largest  value  of  AH^(0  ),  obtained  for  MgO, 

presumably  the  most  ionic  of  the  crystals  treated,  a  value  of  202,3 

kcal/mole  is  obtained.  These  values  are  believed  to  be  more  accurate 

than  earlier  values  given  by  Morris  ( 19 and  by  Huggins  and  Sakamoto  (1957) 

who  obtained  210  ♦  6  and  221  ±  15  kcal/mole.  The  anomalously 

2_ 

low  value  calculated  for  AH®(0  )  for  ZnO  is  believed  to  result 

from  a  substantial  covalent  contribution  in  the  Zn-0  bond  in  this 
oxide. 


IV.  Recommendations 

The  results  of  our  analysis  indicate  that  further  or  continued  work 
should  be  carried  out  in  several  areas  relating  to  equations  of  state 
of  rocks  and  minerals.  These  research  areas  are  outlined  below: 

A.  Strength  effects  at  low  and  high  pressures.  Sew  dynamic  data 
describing  the  stress  differences  supported  by  earth  materials 
upon  passage  of  intense  stress  waves  should  be  obtained.  Knowledge 
of  the  stress  differences  supportable  by  a  j»edium  is  of  Importance 
for  the  following  two  reasons:  I)  The  value  of  the  difference 
In  principal  stresses  present  under  dynamic  conditions  are  critical 
to  the  reduction  of  Hugonlot  or  release  adiabat  shock  data  and  to 
formulation  a  complete  hydrostatic  equarion  of  state  of  the 
material,  2)  Stress  differences  supported  bv  the  media  under  dynamic 
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conditions  directly  enter  into  the  numerical  calculations  of  the 
stress  wave  propagations  and  are  needed  to  formulate  a  complete 
rheological  model  of  the  material  under  varying  dynamic  conditions. 

Further  work  along  the  lines  suggested  by  Walsh  and  Brace 
(1969)  and  Cherry  (1968)  involving  correlation  of  data  from  static 
and  quasi-static  strength  tests  with  dynamic  data  would  aopear 
profitable.  We  have  found  that  the  comparison  of  a  hydrostatic 
equation  of  state  inferred  from  high-pressure  and  high-temperature 
ultrasonic  data  are  also  useful  in  evaluating  stress  differences 
at  high  pressures  (Davies  and  Ahrens,  1969).  Comparison  of  piston- 
cylinder  measurements  with  shock-wave  data,  such  as  is  performed 
by  Peterson  et^  al  (1969)  appears  to  be  another  promising  technique 
to  evaluate  strength  effects. 

Another  approach  to  the  strength  problem  is  the  performance 
of  small-scale  one-dimensional  intense  stress  wave  propagation 
experiments  in  spherical  or  planar  geometry.  In  these  experiments 
one  or  more  of  the  critical  flow  of  variables  are  measured  as 
a  function  of  time  and  thus,  for  example,  pressure  or  particle 
velocity  profiles  are  obtained.  The  various  rheological  models 
for  a  real  material  may  then  be  tested  by  comparing  the  observed 
flow  variable  profiles  with  numerical  calculations.  Some 
experiments  along  these  lines  have  been  carried  out  by  C.  Godfrey 
and  coworkers  (personal  communication,  1968)  in  spherical  geometry. 

Another  important  area  in  which  experimental  data  are 
especially  lacking  is  that  of  the  effect  of  water  content  and 
porosity  on  dynamic  yielding  and  the  stress  differences  supportable 
under  dynamic  conditions. 
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B.  Equation  of  State  Data 

Although  knowledge  of  the  principal  Hugoniots  of  rock  forming 
minerals  is  far  from  complete  and  should  be  continued,  a  considerable 
body  of  data  is  now  available  for  the  oxides  and  silicates.  Data 
pertaining  to  testing  various  mixing  laws  for  employing  equation 
of  state  data  for  minerals  and  oxides  to  predict  the  response 
of  rocks  of  arbitrary  composition  is  important.  Some  work  along 
these  lines  has  recently  been  carried  out  by  McQueen  (1968)  and 
co-workers.  The  effect  of  water  or  a  concentration  of  heavily 
hydrated  minerals  on  the  equation  of  state  should  also  be  examined 
experimentally  and  theoretically  in  future  work. 

C.  Gruneisen's  Ratio 

The  availability  of  thermodynamic  and  ultrasonic  data  at 
high  temperatures  and  moderate  pressures  for  individual  mineral 
species  is  important  in  fixing  parameters  such  as  Gruneisen's 
ratio  under  standard  conditions.  As  recently  pointed  out  by 
several  authors  including  Anderson  (1968)  and  Knopoff  and  Shapiro 
(1959),  knowledge  of  the  Gruneisen  parameter  of  high  pressure  is 
vital  to  the  calculation  of  states  lying  away  from  the  principal 
Hugoniot  as  well  as  a  formulation  of  a  complete  equation  of  state. 
Judging  from  results  of  our  analysis  of  the  high-pressure  data 
for  Gruneisen's  ratio  data  for  stishovite  (Ahrens  et_  al^,  1969a) 
shock  wave  experiments  on  porous  minerals  such  as  those  of  Jones 
et  al  (1968)  on  Sit^  and  Carter  et  al  (1969)  on  MgO  appear  to  be 


a  promising  source  of  data  for  Gruneisen's  ratio  at  high  pressure. 

In  this  regard  the  shock  wave  measurements  themselves,  as  well  as 
the  preparation  of  uniform  samples  and  the  precise  measurements  of 
their  initial  density,  present  challenging  experimental  problems. 

The  new  techniques  in  which  a  large^but  known^ pulse  of  energy 
is  radiatively  deposited  in  a  sample  and  the  rise  in  pressure  is 
then  measured  may  also  be  a  useful  tool  for  generating  new  Gruneisen 
ratio  data  (Graham  and  Hutchison,  1967). 
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Table  3 

Comparison  of  release  adiabat  zero-pressure  and  ultrasonic 

measurements  of  bulk  modulus  K  and  (8K  /3P)  ,  P  ■  0  for 

s  s  s 

various  values  of  initial  strength  correction  La  and  pressure 
coefficient  a  for  MgO. 


T 

(°K) 

Ultrasonic^ 

La  - 

0 

a  ■ 

15  kb 

0 

La 

0 

a 

-  15  kb 

-  0.04 

La 

o 

a 

-  28  kb 

■  o_ 

K 

s 

(Mb) 

1  s  1 
\3P 

S 

K 

s 

(Mb) 

GW. 

K 

s 

(Mb) 

T5T 

s 

K 

s 

(Mb) 

GW. 

300 

1.628 

3.85 

1.693 

3.913 

1.626 

3.913 

1.620 

4.056 

400 

1.612 

3.83 

1.677 

3.926 

1.612 

3.922 

1.608 

4.059 

500 

1.595 

3.81 

1.659 

3.939 

1.596 

3.931 

1.595 

4.063 

600 

1.576 

3.78 

1.640 

3.952 

1.579 

3.941 

1.580 

4.067 

700 

1.556 

3.76 

1.619 

3.966 

1.562 

3.950 

1.564 

4.071 

800 

1.536 

3.74 

1.600 

3.979 

1.544 

3.959 

1.549 

4.074 

900 

1.516 

3.72 

1.581 

3.988 

1.529 

3.964 

1.536 

4.073 

1000 

1.496 

3.71 

1.564 

3.993 

1.514 

3.966 

1.524 

4.070 

1100 

1.476 

3.70 

1.549 

3.998 

1.512 

4.066 

1200 

1.456 

3.69 

Table  4 


Comparison  of  the  pressure  dependence  of  shear  strength 
obtained  from  shock  wave  analysis  ,  with  that  Inferred  from 
ultrasonic  data  according  to  the  hypothesis  of  Towle  and 
Riecker  (1968). 


Substance 

e  (shock) 
s 

e  (ultrasonic) 
u 

MgO  (crystal) 

2.0 

2.0 

A120^  (crystal) 

1.9 

1.1 

(ceramic) 

0.4 

1.1 

Table  5 

Equation  of  State  Parameter 
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FIGURE  CAPTIONS 


Figure  1. 


Figure  2. 


Figure  3. 
Figure  4. 


Pressure-density  relations  for  material  undergoing  a  shock- 

induced  phase  change.  Low  pressure  regime  represents 

£ 

Hugoniot  of  material  with  zero-pressure  density,  Pq. 

High  pressure  regime  represents  Hugoniot  of  material  which 
has  been  completely  converted  to  high  pressure  phase  with 
zero-pressure  density,  pQ.  Metastable  Hugoniot,  adiabat, 
and  isotherm  for  high  pressure  phase  are  centered  at 
state  P  -  0,  T  -  25°C. 

Pressure-density  Hugoniot  curves  for  solid  and  porous 
(sandstone)  a-quartz  and  fused  quartz  in  the  stishovite 
regime.  Calculated  adiabat  centered  at  standard  conditions 
and  25°C  isotherm  for  stishovite  are  calculatd  from  solid 
a-quartz  data.  Isothermal  X-ray  data  for  stishovite 
compares  closely  with  calculated  25°C  isotherm. 

Gruneisen  ratio  versus  density  for  stishovite. 

Lines  of  constant  thermodynamic  gamma  in  the  pressure- 
volume  plane.  Adiabats,  calculated  from  the  data  of 
Spetzler  (1969b),  centered  at  ambient  pressure,  and 
temperature  ranging  from  300°K  to  1200°K  are  indicated. 
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THERMODYNAMIC  VALUE 

SANDSTONE  (JONES  etal,l968)-a 
QUARTZ 

SANDSTONE -METHOD  TL 

FUSED  QUARTZ  (SHIPMAN,  1969) -a 
QUARTZ 

FUSED  QUARTZ  -  a  QUARTZ 
(McQUEEN,  1968) 

FUSED  QUARTZ  (WACKERLE,  1962 
AND  SHIPMAN,  1969)- SANDSTONE 
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